Survival-Time Distribution for Inelastic Collapse 
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In a recent publication [PRL 81, 1142 (1998)] it was argued that a randomly forced particle 
which collides inelastically with a boundary can undergo inelastic collapse and come to rest in a 
finite time. Here we discuss the survival probability for the inelastic collapse transition. It is found 
that the collapse-time distribution behaves asymptotically as a power-law in time, and that the 
exponent governing this decay is non-universal. An approximate calculation of the collapse-time 
exponent confirms this behaviour and shows how inelastic collapse can be viewed as a generalised 
persistence phenomenon. 



There is currently considerable interest in the statisti- 
cal properties of non-equilibrium systems. At the meso- 
scopic scale, the dynamics can often be described by 
Langevin equations and much information can be 
gained from a study of the first-passage statistics of the 
underlying stochastic process 0. Quantitative measures 
such as the survival time or the persistence probability 
H have been introduced to characterise the resistance 
to temporal fluctuations. In many cases the persistence 
probability is found to decay as a power-law in time, and 
the persistence exponent has been determined, using ei- 
ther exact or approximate analytic techniques, for a vari- 
ety of systems. These include the diffusion equation with 
random initial conditions [B, reaction-diffusion systems 
0, phase-ordering kinetics |(| and interfacial growth [Q. 
The persistence exponent has also been measured experi- 
mentally for breath figure s B , a 2-d liquid crystal system 
P], and a 2-d soap froth [ ^0| . 

In a recent publication pd| it was shown that a ran- 
domly forced particle which collides inelastically with a 
boundary can undergo a collapse transition. Namely, if 
the coefficient of restitution is small enough the particle 
will collide an infinite number of times and come to rest at 
the boundary in a finite time. This transition represents 
a novel aggregation mechanism in a driven dissipative 
system and has applications ranging from Brownian mo- 
tion in colloids Jljj to driven granular media |H| . In this 
note we will address the following question: given that 
the particle will collapse at some finite time t, what is 
the distribution of collapse times P c (t), or, alternatively, 
what is the probability that the particle will survive up 
to time t. Surprisingly, we find that the collapse-time 
distribution has power-law tails and that the exponent 
governing the asymptotic decay depends continuously on 
the coefficient of restitution. Furthermore, we are able 
to make a connection between the collapse transition and 
the problem of persistence with partial survival |14j . Our 
results show that inelastic collapse can be viewed as a 
generalised persistence phenomenon. 

We will first outline the arguments which can be used 
to predict the collapse transition, as they will form the 
basis of a systematic calculation of the collapse-time dis- 



tribution. Consider a particle moving in one dimension 
and subject to a random force. The equation of motion 
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where T)(t) is Gaussian white noise with correlator 
{r](t)r)(t')) = 2DS(t - t'). Whenever it returns to the 
origin with speed Vi it is reflected inelastically with a 
reduced speed vj = rv- L , r being the coefficient of resti- 
tution (0 < r < 1). The statistics of the motion of the 
inelastic particle can be inferred from an exact mapping 
onto an elastic problem, which is statistically equivalent 
to a free, randomly accelerated particle. This can be 
achieved by noting that the equation of motion, Eq. ([!]), 
is invariant under the following rescaling of variables 



/ —3 
x = r x; 
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t' = r~H, 



(2) 



whilst in the primed coordinates the velocity is increased 
by a factor 1/r, v' = r~ x v. Thus, the combination of 
an inelastic collision followed by a rescaling of variables 
results in an elastic collision in the new coordinates. If 
one now performs the rescaling after every collision, the 
elapsed time in the inelastic problem can be expressed as 
an integral over the effective elastic variables 



t 



2n(s) 



ds , 



(3) 



where t is the elapsed elastic time and n(s) is the number 
of zero crossings up to time s of a free, randomly acceler- 
ated particle. It was shown in [jll| that if r < r c — e -71 '/^ 
the integral in Eq.(|j|) will converge with probability one 
as t — ► oo. This is the signature of inelastic collapse, as 
the particle will collide an infinite number of times and 
come to rest at time t. Here we are interested in the 
distribution of t (for r < r c ) which is generated by the 
fluctuations in n(t). 

In Figure 1 we plot the collapse-time distribution for 
different values of r, as determined from numerical simu- 
lations in which we evolve many trajectories and form a 
histogram of the resulting collapse times. The details of 
the simulation technique for the inelastic dynamics are 
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given in It is convenient to work with the variable 
T = lni, and Figure 1 shows that asymptotically, the 
probability for collapse within the interval [T, T + AT] 
decays exponentially with a rate which depends on r. 
Thus, the collapse-time distribution, P c (t), does indeed 
show power-law scaling of the form 

P ^)~^jry t- 00 . ( 4 ) 

which serves to define the collapse time exponent 9 c {r). 
Consequently, the survival probability, which is the prob- 
ability that the particle has not collapsed up to time t, 
will decay asymptotically as t~ 8c ( r \ 

One can argue for the value of 9 C in two limiting cases. 
For r = 0, the particle will collapse at its first return to 
the origin because the reflected velocity will be zero. As 
the motion up to that point is just that of a free, ran- 
domly accelerated particle, 9 C (0) can be determined from 
the first passage distribution of the random acceleration 
process. It is known [n5| that this distribution decays 
asymptotically as l/i 5 ^which implies 9 C {0) = 1/4. On 
the other hand, for r > r c there is no collapse transition, 
so that 9 c (r c ) = 0. As can be seen in Figure 1, 9 c (r) 
varies with r between these two limiting values. 
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FIG. 1. The collapse probability P c plotted on a semi- log 
scale as a function of T. Five values of r are shown: 
r = 0.001, 0.01, 0.04, 0.08, 0.12 in ascending order on the right 
of the figure. 

We now turn to an analytic calculation of 9 c (r) which, 
although approximate, does reproduce the simulated val- 
ues rather accurately. The starting point is the integral 
expression for the collapse-time, Eq.(||), where we have 
assumed r < r c and take the limit t — > oo. The lower 
limit of the integral and the initial condition of the par- 
ticle are arbitrary as we are primarily interested in the 
long-time behaviour. For calculational clarity it is con- 
venient to set to = 1 in all that follows jig ]. As n(t) only 



changes at the zero crossings of the underlying elastic 
process, the integral, Eq.(|3j), can be expanded as 

t = h + r 2 (t 2 - ti) + r 4 (i 3 - i 2 ) + . . . , (5) 

where U is the time of the i th crossing. For the random 
acceleration process, the mean number of crossings up to 
time t grows logarithmically with t. Therefor, one can 
map the problem onto a stationary Gaussian process in 
the new time variable T = lnt Jul. It is also useful to 
express the new crossing times Ti in terms of the inter- 
vals between crossings, Tj, of the stationary process. One 
then writes T± = tx, T% = T\ + r 2 . . . and so on. In- 
serting these changes of variable into Eq.(||), and after 
rearrangement one finds 

t = e Tl [1 - r 2 + rV 2 [1 - r 2 + rV 3 [1 + . . . (6) 

Next we make an approximation which has been em- 
ployed in similar first-passage and persistence problems 
P,p7{ , the independent interval approximation (IIA) . We 
assume that the are all independent and drawn from 
the same distribution Pj(t), the interval size distribution. 
Within this approximation, and exploiting the infinite 
hierarchical structure of Eq.(^), the stochastic variable t 
can be expressed as 

t = e T [l-r 2 + rH'], (7) 

where t' is a random variable drawn from the same dis- 
tribution as t, P c (t). It follows that P c (t) satisfies the 
integral equation 

/OO /"OC 
dfP c (t') / dTPi(T)S(t-e T [l-r 2 +rH']). 

(8) 

We are interested in the asymptotic behaviour of P c (t) 
which we expect to have the form P c (t) ~ l/t 1+ec . In- 
serting this into Eq.(||) and performing the asymptotic 
analysis, we are left with the following equation which 9 C 
must satisfy, 

y* oo 

r 2e - / dTe 6 ° r Pj(t) = 1. (9) 
Jo 

Thus, from knowledge of the interval size distribution 
Pi (r) one can determine the collapse-time exponent 9 C . 

The random acceleration process has been studied by 
Burkhardt |t5|, and the interval size distribution has 
been calculated up to quadrature for arbitrary initial 
conditions. Asymptotically it is known to decay expo- 
nentially, having the form P/(r) ~ aie~ r / 4 . The coef- 
ficient ai can be determined exactly from the results in 
p5| . The survival probability, up to time t, for a particle 
injected from the origin with speed uq has the asymptotic 
form 

Q(0,u ,t)= 3 |u f . (10) 
V27rr(3/4) ti 
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This must be averaged over uq drawn from the scaling 
distribution and correctly weighted to pick out the ze- 
ros of x; the resulting distribution for u is P(uq) = 
2(itoAo)e~ u °/ to , where to is the injection time and P(uq) 
is normalised over the interval [0, oo]. Performing the 
average over uq in Eq.fJToJ), one finds that the interval 
size distribution has the expected form aie~ T / 4 , where 
t = hit /to and 



a i 



3 T(5/4) 
4 v / 2^r(3/4) 



= 0.221314. 



(11) 



Figure 2 shows Pi(t) determined numerically and con- 
firms the exponential decay at large r with the correct 
value of ai given above. 
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FIG. 2. The interval size distribution Pi(t) plotted on a 
semi-log scale as a function of r. The points are the numeri- 
cal data while the line shows the approximate fit, Eq.(12). 

Knowledge of Pi(t) for all r is needed to evaluate the 
integral in Eq.(^). A good fit to the numerical data can 
be obtained by writing Pi(t) as the sum of just two ex- 
ponentials, 



a 2 e 



(12) 



with 9i = 1/4 and a% = 0.2213. The parameters 0i and 
<i2 can be determined by imposing two constraints: (1) 
that Pi(t) is normalized over the interval [0, oo] and (2) 
that it predicts the correct density of zeros U which is 
governed by the first moment of Pj(r), p = V3/2tt — 
1/(t). These conditions fix the unknown parameters to 
be 02 = 1.3254 and 02 = 0.1521, and this approximate 
expression is also plotted in Figure 2. It can be seen that 
the parameterised -Pj(t) fits the simulation data rather 
well over the whole range of r. 

When this analytic expression for Pj(t) is inserted into 
Eq.(||) one finds the following transcendental equation, 
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a 1 



a-2 



(13) 



which can be used to determine 6 c (r). The numerical so- 
lution of this equation is shown in Figure 3, together with 
the exponent values measured directly from the simula- 
tions. For r close to r c it is hard to access the asymptotic 
regime numerically, hence the lack of data in this region 
of the plot. However, there is rather good agreement for 
all the values of r considered. Note that with our param- 
eterised expression for Pi(t), the IIA correctly predicts 
both the r — * limit and the critical value of r, namely 
r c = e -71 "/^, where 9 C — > 0. Thus, the analytic calcula- 
tion confirms that the collapse exponent is non-universal 
and depends continuously on r. The inset in Figure 3 
shows the behaviour of 6 c {r) for r — ► 0, where the ap- 
proximate analytic calculation and the simulation data 
appear to coincide. To leading order in r, Eq.(|l^) can be 
expanded to give 



1/4- aiV7 + 



(14) 



As only the first term on the right hand side of Eq. (|13j ) 
contributes in this limit, the expansion is exact within 
the IIA, since both #i(= 1/4) and a\ are known exactly. 
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FIG. 3. The collapse-time exponent 9 C as a function of r. 
The points are the numerical data, obtained from measuring 
the slopes of P C (T). The errors are of the same order as the 
size of the circles. The line is the analytic approximation de- 
termined from the solution of Eq.(13). The inset shows an 
expanded view of the main figure for r close to zero where 
the simulation data and the approximate calculation appear 
to coincide. 

Finally we point out an intriguing connection between 
the exponent 9 C and the persistence exponent of a ran- 
domly accelerated particle with 'partial survival' Jl4) . 
Consider the process x(t) described by Eq.(|l]), and let the 
particle survive with probability p each time its position, 
x(t), crosses through zero. Working in the logarithmic 
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time variable, T = hit, the probability that the particle 
survives for a time T (in the stationary state) is 



Psurv(T) 



n=0 



p n Pn(T), 



(15) 



where P n (T) is the probability that the process has 
n zero crossings in time T. We anticipate |t4| that 
P surv (T) ~ e - e ^ T for T -> oo, i.e. that the Laplace 
transform P S urv(s) will have a simple pole at s = —9{p). 

The Laplace transforms P n (s) were worked out, in 
terms of the Laplace transform P/(s) of the interval-size 
distribution, Pj(r), within the IIA in [Q. Carrying out 
the sum over n in Eq. (|l5|) gives -P SU rv(s), and one finds 
a simple pole given by the condition pPi(s) — 1. Setting 
s = -0(p) yields finally p J™ dr e e ^ r Pj(t) = 1 as the 
equation determining This equation has an iden- 

tical form to Eq.(||) for 9 C , suggesting that the collapse 
process with coefficient of restitution r is in some sense 
equivalent to a randomly forced particle with survival 
probability p — r 20c at each collision. As yet, however, 
we have been unable to establish this equivalence outside 
the IIA. 

To summarize, we have studied the collapse-time dis- 
tribution, P c (t), for a randomly- forced particle colliding 
inelastically with a boundary, in the regime r < r c where 
inelastic collapse occurs. We have shown that the distri- 
bution has a power-law tail, P c (t) ~ \/t 1+0c , character- 
ized by an exponent 9 C which varies continuously with r 
for < r < r c . An analytical approach, based on the ap- 
proximation that the intervals (in the variable T = hit) 
between zero crossings of a free, randomly accelerated 
particle can be treated as statistically independent, leads 
to results which are in good agreement with numerical 
simulations. We have also noted a formal similarity be- 
tween Eq.(^), which determines 9 C within this approach, 
and the equivalent equation for the persistence exponent 
of a randomly accelerated particle with partial survival 
ijLH . This similarity invites further study. 

This work was supported by EPSRC (UK) under grant 
number GR/K53208. 
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